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ABSTRACT 

Comprehensive vibration-rotation line lists for eight isotopologues of carbon mono- 
sulphide (CS) ( 12 C 32 S, 12 C 33 S, 12 C 34 S, 12 C 36 S, 13 C 32 S, 13 C 33 S, 13 C 34 S, 13 C 36 S) in 
their ground electronic states are calculated. These line lists are suitable for temper¬ 
atures up to 3000 K. A spectroscopically-determined potential energy curve (PEC) 
and dipole moment curve (DMC) are taken from literature. This PEC is adapted to 
suit our method prior to the computation of ro-vibrational energies. The calculated 
energies are then substituted by experimental energies, where available, to improve the 
accuracy of the line lists. The ab initio DMC is used without refinement to generate 
Einstein A coefficients. Full line lists of vibration-rotation transitions and partition 
functions are made available in an electronic form as supporting information to this 
paper and at www. exomol. com. 

molecular data; opacity; astronomical data bases: miscellaneous; planets and satellites: atmospheres; stars: low-mass 


1 INTRODUCTION 


It was the development of radio astronomy that led to the realisation that the majority of our galaxy, and hence the universe, 
is dominated by molecular processes. More than 150 molecules have been detected so far in t he interstellar medi um (ISM) 
by direct observation of their spectra. Carbon monosulphide (CS) is one of these molecules (IPenzias et al.l 119711 1 and is a 
diatomic of both atmospheric as well as astrophysical interest. In th e Earth’s atmo sphere, CS plays a role in the formation of 
aerosols, in particular carbonyl sulphide (OCS) in the tropos phere dLi et alj 201^1. 


In t he solar system C S has been observed in comets IlCanaves et al 


20071 and the collision of comet Shoemaker- 


Levy 9 I Orton et al. 199511 led to its detection in the atmosphere of Jupiter jMoreno et al.l 20031). Astronomically the 
molecule has been observed in a variety of objects suc h as carbon-rich stars i Breerman et al . 197Sl : [Botsdrwim^^Jgebaldl 


I Henkel & Ball^ 1985^ 


1985^ Agundez&Cemicharo 2006|h star f orming regions ( Davis et, al. 20131 ) and dense interstellar clouds |Niisso7e7aLl 200(1 : 
McQuin n et al. !7xh Scoyilleet alJ 201 51). In fact CS is o ne of the most abundant sulphur-containing species in interstel- 


lar clouds dshi_et__al, 2011 : iBilalbegovic fe Baranovic 2015 ) with sever al isotopologues long detected outside the Milky Way 


Maucrsbcrger et al.ll 1989b 


aj; 


Henkel et al. 


19931) . 


The numerous astronomical detections of CS and the importance of the molecule in our own atmosphere has motivated 
copious laboratory studies. Experimentally the CS spectrum has been st udied in wavelengt h r egions ranging from the mi¬ 
crowave to the ultraviolet (UV). A pioneering study was carried out by 1 Crawford fe Shurclifli (1934)) whom discovered the 


main A 1 II - X 4 S + transit i on in t he visible. Add i tional elect ronic transit i ons in the visible to near UV have been investigated 
by, for example, 


Bell et al. 


( 19721) , Cossart et al. ( 1977 ) and Stark et al. ( 1987 ). 


Mockler fe Bird ( 1955 ) made the first measurements of rotational lines for the 12 C 32 S, 12 C 33 S, 12 C 34 S and 13 C 32 S iso - 


topomers in the microwave region. Transitions in this region have also been observed by, for example, Lovas fe Krupenlel (Il974 ). 
Early work in the millim etre wave region b egan with measurements of rotatio nal 12 C 32 S and 12 C 34 S lines by iKewlev et al. 


1963 ). later extended by Bogey et ali ( 19821 ) who additionally observed 13 C 32 S. Bogey et al. ( 1981 ) also presented a study of 


the millimetre spectrum of rarer isotopologu es 12 C 33 S, 12 C 36 S, 13 C 33 S and 13 C 34 S. More recent work in t he region has been 


carrie d out bv lAhrens fe Winnewisse? [|l990ll f 12 C 32 S, 12 C 34 S, 12 C 33 S, 13 C 32 S, 12 C 36 S, 13 C 33 S, 13 C 34 S). Eim fc Yamamotcl 
( 20031) f 12 C 32 S. 12 C 34 S) and lGottlieb et al.l I 2 OO 3 I) ( 12 C 32 S, 12 C 34 S, 12 C 33 S, 13 C 32 S). 


The first study in the infra-red region was performed by Todd (19771 ) who measured the v = 2 — 0 vibrational band of the 
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mam 
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isotope while Todd Sz Olson ( 19791 ') and Yamada fc Hirota ( 19791 ) measured several Aw = 1 bands of 12 C 32 S, 12 C 33 S, 
C 34 S and 13 C 32 S. The r esulting molecular parameters were later refined by Winkel et al. ( 1984) and Brujcholder_et_alJ 
1 1987 ). Winkel et al. 1 1984 ) measured many Av = 2 bands for vibrational levels up to v = 8 while Burkholder et alT ( 1987 ) 
obtained high resolution measurements of the 1-0 band, and 2-1 band for the main isotope, of 12 C 33 S, 12 C 34 S and 13 C 32 S for 
J up to 41, 28, 32 and 28 respectively. Aw = 1 bands up to v = 9-8 for the main isotope were l a ter m easured by Ram et ahl 

who reported Aw = 
7-6 to a higher 


<ioTi), 


(119951) . The most recent research on CS infra-red spectra has been performed bv lUehara et al 
1 transitions of 13 C 32 S up to w = 5-4. Additionally they measured Aw = 1 transitions of 12 C 32 S up to v 


up to v = 
1. (1 19951) . 


accuracy than lRam et al. 

Using the vibration-rotation and pure rotation data on this molecule available to them, Coxon fc Haiigeorgiou ( 19921 ) 
derived a spectroscopic potential energy curve (PEC) that reproduced all the input experimental data within experimental 

error. This PEC is the starting point for the current work. _ _ _ 

Transition probabilities or Einstein A coefficients have been provided by Botschwina fe Sebald ( 19851) and Pineir o et al. 


( 19871) . The latter produced transition lists for rotational quantum numbers J' — J = ±1 for J ^ 200, and vibrational w ^ 20, 
though only for the four most abundant isotopes ( 12 C 32 S, 12 C 33 S, 12 C 34 S and 13 C 32 S). Line lists for all isotopologues of CS, 
including some rotation-vibrati on transition s, are a vailable from CDMS jMiillet_^t_^iL 2005 ). These we r e constructed from 
experimental data obtained by Bogey et al. |l982 h_ Boge^j^al. J 1981 ). Ahrens^fc^Vimewissei 1 19991) . Kim fc Yarriamoul 


20031) , Gottlieb et al. 1 20031) . Burkholder et al.l~ 1987). 


Ram et al. 


Winkel et al. 


and the only available exper¬ 
imental measurements of the dipole moment (IWinnewisser fc Cook! 119681 ). The line lists are very accurate and recommended 
for use in radio astronomy, though are limited to w < 4 and J ^ 99 and hence temperatures, T below about 500 K. The aim 
of this work is to produce comprehensive line lists for all stable isotopologues of CS suitable for modelling hot environments, 
such as carbon stars (T ~ 3000 K). 

The ExoMol project aims to provide line lists on all the molecular transitions of importance in the atmospheres of planets 
I Tennyson fc Yurchenko 2Q12 ). T he Ex oMol methology ha s already been a pplied to a numbe r of di atomic molecules: BeH, 
MgH and CaH (I Y adiri e t al. 20121), Si O ([Barton et al 'l2013l) . NaCl and KC1 ( Barton et al. 2014 ). PN ( Yorke et al.ll2014 ). AlO 
( Patrascu et al. 20151) and NaH ( Rivlin et al. 2015 ). In this paper, we present ro-vibrational transition lists and associated 
spectra for all stable isotopologues of CS. 


2 METHOD 


The line lists for all eight isotopologues of CS, which we have named JnK, were obtained by solvin g the Sch rodinger equation 
allowing for Born-Oppenheimer Breakdown (BOB) effec ts using the program LEVEL 8.0 ( Le Roy 200?! ). In principle the 
calculations were initiated using the spectroscopic PEC of Coxon fc Haiigeorgiou ( 19921 ). In practice, as described below, the 
PEC was first expressed in a form compatible with LEVEL. The PEC was then adapted to improve the results of calculations 
performed using LEVEL. After computation the line lists are improved by replacing calculated energies with experimentally 
derived energies, where available, and shifting th e remaining energies to maintain LEVEL predicted energy level separations. 
A theoretical dipole moment curve (DMC) from I Pineiro et alj (119871) was also employed. 


2.1 Potential Energy Curve 


We did not g enerate a new PEC for CS. A full set of potential parameters representing a very accurate PEC is already 
available from Coxon fc Haiigeorgiou 1 19921) . The authors employed data for four isotopomers ( 12 C 32 S, 12 C 33 S, 12 C 34 S and 
13 C 32 S) in a least-squares fit to a potential function in the Born-Oppenheimer approximation and BOB functions to determine 
a PEC valid for all isotopologues. They set the dissociation energy D e to 59300.0 cm -1 and expressed their fitted potential 
Uqs ( R ) as a Morse potential with a variable /3: 


uEs(R) = D e (1 - exp [-0(R){R - Re)}) 2 , 
where: 

0{R) = Ho + pl {R — Re) + 02 (i? - Re) 2 + ... 

The isotopically invariant breakdown functions were modelled as: 

U c {R) = T, i=1 u < t{R- R e )\ 

and 

U s {R) = E i=1 uUR~ Re)\ 

while the J-dependent non-adiatbatic breakdown function was modelled as: 


(1) 


( 2 ) 


( 3 ) 


( 4 ) 
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qcs(R) = q c (R)/M c + q s (R)/M s , 


(5) 


where: 

q (k \R) = ^ i= ipi k) (R-R e y, (6) 

such that q( k \R e ) = 0. These were applied in the effective radial Hamiltonian according to: 

Hcs(R) = + U& f (R) + ( Plt/R 2 )J{J + 1)[1 + qcs(R)}, (7) 

where /3 2 t = ti 2 /2p a t, defined with the atomic masses and: 

ug (R) = uEs (R) + Uc (R)/Mc + U S (R)/Ms , (8) 

The variable f) Morse is not implemented in LEVEL8.0. The coefficients from Coxon fe Haiigeorgiou ( 19921 ). given in Column 
II of Table [TJ were hence used to generate turning points for a range of internuclear distances. The turning points were used 
directly in LEVEL. 

Employing the potential parameters of Coxon fe Haiigeorgioul (1 199211 in thi s way we could not reproduce the vibrational 
energies to the spectroscopic accuracy achieved by Coxon fe Haiigeorgioul ( 19921 ). By applying small ‘corrections’ to potential 
parameters uf, uf, pf and p®, we were able to predict the ro-vibrational energies up to v = 9, and experimental frequencies, 
for 12 C 32 S and 13 C 32 S to within 0.02 cm -1 and 0.04 cm” 1 respectively (see JnK columns in Tableland Tabled. To give an 
example, the residual (obs-calc) for 12 C 32 S T\ without the ‘corrections’ was 0.027 cm” 1 and 0.007 cm” 1 after correction. The 
term ‘correction’ is used tentatively in this context as, although the modificati on of the potential parameter s improved the 
present results, this is not an improvement on the variable (3 Morse presented by Coxon fc Haiigeorgiou ( 1992l h The potential 
parameters used in this work are given as Column III of Table [T| 

To further improve our results we took advantage of the ExoMol format used to store the line list, see Section 3. Put simply 
this is a states file containing level energies and a transitions file detailing allowed energy level couplings. The advantage of the 
format is it gives the option of replacing calculated energies with more refined or experime ntal energies s uch th at, when the 
files are unpacked to produce the line list, more accurate line frequencies are computed, see Barber et al. ( 20141 ) for example. 

First we att empted to refine our ro-vibrat ional energies (E 3 "^) for all isotopologues using the vibrational (J = 0, v ^20) 
energies given in [Coxon fe Haiigeorgioul (ll992l ) ( E ^q x ) and the formula: 


e: 


JnK—Cox 
v,J 


_ / TT'JnK rpJnKX I T7>CoX 

— \^v,J — &v,0 ) ~T EU v ,0 


(9) 


Ro-vibrational energies for v > 20 were shifted to maintain the energy level separations predicted by LEVEL according to: 


TiJnK—Cox 
v,J 


_ / T-iJnK T-iJnK \ . 7 —iCox 


( 10 ) 


We co uld then reproduce the vibrational energies to the same spectroscopic accuracy achieved by Coxon fe Haiigeorgiou 
1 19921b nowever not all the line frequency predictions improve d (see JnK-Cox column s in Ta ble Q). 

This is likely due to the fact experimental da ta available to Coxon Ik Haiigeorgiou ( 1992l l was limited to J ^41 and J $)28 
while Ram et al. (1995) and Uehara et al. 1 2015h assigned lines for J up to 113 and 86 for 12 C 32 S and 13 C 32 S respectively. 
Therefore we d ecided to determine experimental energies directly from frequencies measured by Ram et al. d 1995h and 


Uehara et al. 1 2015h using the measured active rotation-vibration energy levels (MARVE L) technique ( Furtenb acher et al 
2007 I 1 which involves inverting transition frequencies to extract experimental level energies. lUehara et al l J2015 b is the more 


accur ate experimental s tudy and thence energies extracted from these frequencies were used preferentially over those extracted 
from Ram et al. 1 1995li frequencies where possible. We extracted 733 energies in total for the main isotopologue and 341 
energies for 13 C 32 S, see Table [6] 

12 C 32 S and 13 C 32 S energies for the experimental ranges were replaced with experimentally derived energies. Ro-vibrational 
energies for ( v , J) outside the experimental ranges were shifted to maintain the energy level separations predicted by LEVEL 
according to the following equations. For v < an d J > J® 

rpJnK — Exp / 7-iJnK T-iJnK\ , rpExp /1 1 \ 

E v ,j = [Ev.J ~ Ev,o ) + E v 0 ( 11 ) 

For v > « Exp • 

l ^ ^max • 

rpJnK—Exp / rpJnK rpJ n K \ . rpJ n K—Exp /io\ 

E v,j = [Ev,j - Ev-I,j) + E V _ 1}J (12) 

The experimental frequencies for 12 C 32 S and 13 C 32 S, by default, were reproduced almost exactly using this method. The 
largest residual is < 0.001 cm” 1 . Hence the accur acy of these line li sts s hould be equal to the experimental accuracies, which 
are expected to be 0.012 cm” 1 and 0.01 cm” 1 for 


Ram et al. ( 1995 ) and Uehara et al. ( 20151 ) respectively. 

For 12 C 33 S and 12 C 34 S lBurkholder et al . (198711 measured infra-red v = 1 0 absorption frequencies at 0.004 cm” 1 unapodized 


resolution. The Coxon refined energies, see Eqs. I9land ll0l reproduce these frequencies to very high precision (see Table [4j, as 
would be expected considering they fitted to these frequencies. The Coxon refined energies also represent an improvement on 
the unrefined JnK energies (see Tabled, therefore we chose to employ them in our final line lists for these two isotopologues. 
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Table 1. Coefficients of the Born-Oppenheimer potential and radial functions for the X state of CS. 


Parameter 

Coxon & Haiigeorgiou (1992) 1 

This Work 

Re A 

1.5348175 (27) 

1.5348175 

/3o 

1.89830910 (907) 

1.89830910 

01 

0.0173292 (2255) 

0.0173292 

@2 

0.1220304 (5245) 

0.1220304 

03 

0.0317616 (6808) 

0.0317616 

01 

0.0437801 (7351) 

0.0437801 

05 

0.028153 (1245) 

0.028153 

U 1 

-473.55 (11.10) 

-474 

“2 

1631.08 (21.83) 

1631.08 

4 

-1380.3 (640.8) 

-1380.3 

u c 

1 

u s 

-13714.9 (730.8) 

-13714.9 

-448.24 (18.61) 

-448 

u s 

1124.90 (52.07) 

1124.90 


-1245.2 (116.2) 

-1245.2 

4 

P2 

-0.0041274 (3021) 

-0.0038353 

-0.0008036 (4672) 

-0.0003364 


1 Uncertainties are given in parentheses in units of the last digit. 


Table 2. A comparison of theoretically and experimentally derived vibrational term values for 12 C 32 S and 13 C 32 S in cm 1 . 


T v 

Experiment 

Calculated 

Obs-Calc 

Coxon & Haiigeorgiou (1992) 




This Work 

This Work 




(JnK-Exp) 

(JnK) 

(JnK) 

(JnK-Cox) 

(JnK-Cox) 



12 C 32 S 




Uehara et al. (2015) 
Ti 

1272.162085 

1272.1690 

-0.0069 

1272.16214 

-0.00006 

t 2 

2531.353715 

2531.3486 

0.0051 

2531.35384 

-0.00013 

t 3 

3777.597715 

3777.5899 

0.0078 

3777.59800 

-0.00029 

t 4 

5010.916850 

5010.9087 

0.0082 

5010.91739 

-0.00054 

T S 

6231.333760 

6231.3229 

0.0109 

6231.33456 

-0.00080 

n 

7438.870820 

7438.8584 

0.0124 

7438.87183 

-0.00101 

Tv 

8633.550080 

8633.5338 

0.0163 

8633.55125 

-0.00117 

Ram et al. (1995) 

Tg 

9815.39254 

9815.3752 

0.0173 

9815.39457 

-0.0020 

Tg 

10984.42006 

10984.4025 

0.0176 

10984.42297 

-0.0029 



13 C 32 S 




Uehara et al. (2015) 
Ti 

1236.315929 

1236.3524 

-0.0364 

1236.31591 

0.00002 

t 2 

2460.388235 

2460.3782 

0.0100 

2460.39043 

-0.00220 

t 3 

3672.237852 

3672.2234 

0.0145 

3672.24479 

-0.00694 

n 

4871.88567 

4871.8831 

0.0026 

- 

- 

t 5 

6059.35246 

6059.3547 

-0.0021 

- 

- 


For the remaining isotopologues, with th e exception of 13 C 36 S which h as not been observed experimentally, we have only 
measurements of rotational frequencies from Ahrens fe Winnewisser 1 1999li to compare with. These are expected to have an 
accuracy of at least 0.00002 cm” 1 . As can been seen in Table [5] the agreement is excellent. Due to our methods, unrefined 
JnK and Coxon refined energies predict the same rotational frequencies. However, since the Coxon refined energies improved 
ro-vibrational frequency predictions for other isotopologues, these are employed in our final line lists for the remaining four 
isotopologues. 

An overview of the energy level content of our final ’hybrid’ line lists is given in Table [7] Although we have used 
terms JnK, JnK-Cox and JnK-Exp in the text to refer to unrefined, Coxon refined and experimentally substituted energies 
respectively, the final line lists as provided in supplementary data and on www.exomol.com are simply named JnK. 

We note that the behaviour of the potential curve for longer internuclear distances than 3 A was found to b e divergent 
and n on-physical resulting in the m olecule not dissociating prope rly. This is consistent with the findings of I Coxon fc Colin 
11997), that the model employed by Coxon fc Haiigeorgiou (1992) does not account for the inverse-power behaviour of the 
PEC at long range. For this reason we only generated turning-points for internuclear distances between i? m i n = 1.00 A and 
Umax = 3.00 A with a grid spacing of 0.0007 A. This has consequences for the temperature range considered. Based on our 
partition sum, see Section 2.3, this range now extends to 3000 K. 
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Tab le 3. Compa rison of predicted ro-vibrational frequencies, in cm 1 . with experimental line positions measured by I Ram et al.l ll 19951 1 
and lUehara et al.l 1120151) for 12 C 32 S and 13 C 32 S. 


J' 


12 C 32 S 


Experimental 

JnK-Exp 


Calculated 

JnK 


Obs-Calc 

JnK 


Calculated 

JnK-Cox 


Obs-Calc 

JnK-Cox 


Uehara et al. (2015 


10 

9 

1 

0 

10 

11 

1 

0 

9 

10 

2 

1 

11 

10 

2 

1 

100 

101 

2 

1 

20 

21 

6 

5 

21 

20 

6 

5 

78 

79 

6 

5 

10 

11 

7 

6 

15 

14 

7 

6 

63 

64 

7 

6 

102 

103 

1 

0 

107 

106 

1 

0 

107 

106 

2 

1 

89 

88 

6 

5 

89 

88 

7 

6 

25 

26 

8 

7 

30 

31 

8 

7 

52 

53 

8 

7 

59 

58 

8 

7 

25 

26 

9 

8 

28 

27 

9 

8 

52 

53 

9 

8 

59 

13 C 32 S 

58 

9 

8 

2 

1 

1 

0 

79 

80 

1 

0 

12 

13 

2 

1 

70 

71 

2 

1 

7 

8 

3 

2 

53 

52 

3 

2 

6 

7 

4 

3 

64 

65 

4 

3 

14 

13 

5 

4 

45 

46 

5 

4 


1287.847083 

1253.542068 

1242.440398 

1276.248067 

1040.852470 

1172.020927 

1237.819388 

1049.137636 

1176.836843 

1216.680436 

1072.087367 

Ran^et^aL (1995) 

1047.2868 

1371.8550 

1357.5752 

1296.2749 

1282.3246 

1137.7464 

1128.3927 

1084.0478 

1251.2061 

1125.2391 

1207.1842 

1071.8531 

1237.6787 

i Uehani i e^ i aL (2015) 

1239.369011 

1080.972174 

1203.322346 

1089.992782 

1199.383201 

1276.188741 

1188.850787 

1080.169544 

1207.300671 

1107.707636 


1287.8541 

-0.0070 

1253.5490 

-0.0070 

1242.4286 

0.0118 

1276.2363 

0.0118 

1040.8485 

0.0039 

1172.0201 

0.0008 

1237.8182 

0.0012 

1049.1376 

0.0000 

1176.8360 

0.0008 

1216.6786 

0.0018 

1072.0888 

-0.001 

1047.2910 

-0.004 

1371.8529 

0.0021 

1357.5763 

-0.0011 

1296.2749 

0.0000 

1282.3263 

-0.0017 

1137.7446 

0.0018 

1128.3910 

0.0017 

1084.0496 

-0.0018 

1251.2079 

-0.0018 

1125.2367 

0.0024 

1207.1833 

0.0009 

1071.8548 

-0.0017 

1237.6766 

0.0021 


1239.4051 

-0.0360 

1080.9775 

-0.0053 

1203.2753 

0.0470 

1089.9459 

0.0469 

1199.3761 

0.0071 

1276.1977 

-0.0089 

1188.8625 

-0.0117 

1080.1756 

-0.0061 

1207.3057 

-0.0051 

1107.7159 

-0.0083 


1287.8472 

-0.0001 

1253.5421 

- 0.0000 

1242.4407 

-0.0003 

1276.2484 

-0.0003 

1040.8606 

-0.0082 

1172.0219 

-0.0010 

1237.8200 

-0.0006 

1049.1394 

-0.0018 

1176.8400 

-0.0032 

1216.6827 

-0.0022 

1072.0928 

-0.0054 

1047.2841 

0.0027 

1371.8459 

0.0090 

1357.5884 

-0.0132 

1296.2767 

-0.0018 

1282.3303 

-0.0057 

1137.7465 

-0.0001 

1128.3928 

-0.0001 

1084.0514 

-0.0036 

1251.2098 

-0.0037 

1125.2379 

0.0012 

1207.1844 

-0.0002 

1071.8559 

-0.0028 

1237.6778 

0.0009 


1239.3686 

0.0004 

1080.9411 

0.0311 

1203.3239 

-0.0015 

1089.9944 

-0.0016 

1199.3855 

-0.0023 

1276.2070 

-0.0182 

1188.8411 

0.0097 

1080.1542 

0.0154 

1207.3057 

-0.0051 

1107.7159 

-0.0083 


Table 4. Comparison of predicted ro-vibrational frequencies, in cm 1 . with experimental line positions measured bv lBurkholder et alj 
lll987ft for 12 C 33 S and 12 C 34 S. 


J' 

J" 

v' 

v" 

Experimental 

Calculated 

JnK 

Obs-Calc 

JnK 

Calculated 

JnK-Cox 

Obs-Calc 

JnK-Cox 

12 c 33 s 

3 

2 

1 

0 

1271.73537 

1271.736018 

-0.000648 

1271.735366 

0.000004 

4 

5 

1 

0 

1258.72386 

1258.724631 

-0.000771 

1258.723979 

-0.000119 

29 

28 

1 

0 

1308.72699 

1308.728068 

-0.001078 

1308.727416 

-0.000426 

25 

26 

1 

0 

1221.09717 

1221.098153 

-0.000983 

1221.097501 

-0.000331 

12 C 34 S 

3 

2 

1 

0 

1266.78054 

1266.775300 

0.005240 

1266.780744 

-0.000204 

4 

5 

1 

0 

1253.87137 

1253.865751 

0.005619 

1253.871195 

0.000175 

35 

34 

1 

0 

1310.80236 

1310.797693 

0.004667 

1310.803137 

-0.000777 

35 

36 

1 

0 

1197.09619 

1197.091827 

0.004363 

1197.097271 

-0.001081 


© 2015 RAS, MNRAS 000.ITWT21 























6 Paulose et al 


Table 5. Comparison of predicted rotational frequencies, in cm 1 . with experimental line positions measured by I Ahrens &; WinnewisseH 
dl999h for 12 C 36 S, 13 C 33 S and 13 C 34 S. 


J' 

J" 

v' 

v" 

Experimental 

Calculated 

JnK JnK-Cox 

Obs-Calc 

JnK JnK-Cox 

12 C 3B S 

6 

5 

0 

0 

9.507279 

9.507297 

-0.000018 

22 

21 

1 

1 

34.561685 

34.561709 

-0.000024 

13 C 33 S 

6 

5 

0 

0 

9.173973 

9.173971 

0.000002 

20 

19 

0 

0 

30.545844 

30.545846 

-0.000002 

13 c 34 s 

20 

19 

0 

0 

30.293201 

30.293299 

-0.000098 

13 

12 

2 

2 

19.429173 

19.429154 

0.000019 


Table 6. Summary of Energies Extracted from Experimental Frequencies. 


V 

Jmax 

Total Extracted 

Using Uehara et al. (2015) 

Using Ram et al. (1995) 

12 c 32 s 

0 

106 

107 

87 

20 

1 

106 

107 

102 

5 

2 

101 

102 

94 

8 

3 

93 

94 

86 

8 

4 

89 

90 

89 

1 

5 

89 

80 

62 

18 

6 

76 

55 

41 

14 

7 

70 

46 

24 

22 

8 

59 

32 

0 

32 

9 

13 C 32 S 

59 

32 

0 

32 

0 

80 

78 

78 

0 

1 

71 

70 

70 

0 

2 

70 

69 

69 

0 

3 

65 

64 

64 

0 

4 

46 

30 

30 

0 

5 

46 

30 

30 

0 


One approach proposed by Haiigeorgiou fe Le Roy ( 2000l ) to overcome the proble m is to employ a Modified Lennard- 

in|ll 


Jones (MLJ) potential function in place of the Morse variable /3 (|Coxon fc Coli: 
temperature line list without this treatment, it is beyond the scope of this work. 


1997). As it is possible to produce a high 


2.2 Dipole Moment Curve 


The only experimen tal dipole moment data for CS in the literat ure are Stark measurements o f the J = 1 t— 0 matrix elements 
for v = 0 and v = 1 ( Winnewisser fc Cookll l OfiJ . This motivated Botschwina fc Sebaldl ( 19851 ) to perform ab initio calculations 
of the dipole moment for CS over a range of internuclear separations (2.5 ao < R < 3.3ao)- The resulting DMC yielded results 
that were in good a greement with the Stark values and so it was extended to longer and shorter internuclear distances by 


Pineiro et al.l (| 1987 ) using a Pade approximant: 

M(x) = - N ° + NlX - 

1 + D\x + D2X 2 + D3X 3 + -D 43; 4 + D$x 5 


(13) 


Table 7. Overview of sources of energy levels used in the JnK ’hybrid’ line lists. 


Isotopologue 

Experimental Energies 

Coxon Energies 

Shifted Energies 

12 c 32 s 

v ^ 9, J ^ 106 

None 

v > 9, J > 106 using Eg. inland Eg. 1121 

12 c 33 s 

None 

v ^ 7 using Eq. [9] 

v > 7 using Eq. [10] 

12 c 34 s 

None 

v ^ 7 using Eq. |9] 

v > 7 using Eq. [10] 

12 c 36 s 

None 

v ^ 7 using Eq. [9] 

v > 7 using Eq. [10] 

13 c 32 s 

v ^ 5, J ^ 80 

None 

v > 5, J > 80 using Eq. 1111 and Eq. 1121 

13 c 33 s 

None 

v ^ 3 using Eq. [9] 

v > 3 using Eq. [To] 

13 c 34 s 

None 

v ^ 3 using Eq. [9] 

v > 3 using Eq. [10] 

13 c 36 s 

None 

v ^ 3 using Eq. [9] 

v > 3 using Eq. [10] 
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R[A] 

Figure 1. Dipole moment curve for CS obtained from the Pade expansion parameters of lPineiro et al.| i 1987ft . 


Table 8. Comparison of 12 C 32 S partition functions 


T(K) 

This work 

CDMS 

HITRAN 

150 

127.9828 

127.9818 

128.1770 

300 

256.3156 

256.3136 

257.0819 

500 

437.5903 

437.5868 

439.6790 

1000 

1018.57 

- 

1026.38 

2000 

2880.03 

- 

2910.73 

3000 

5705.80 

- 

5803.37 


where x = (R — R e )/R e . This function and coefficients Ni, Di were used to generate turning points for input to LEVEL. The 
results of Pineiro et al. (1 9871 1 were reproduced to the precision quoted in their paper by this approach. 


2.3 Partition Functions 


Partition function values for all 8 isotopologues of CS were calculated by a direct sum of all calculated energies for a range 
of temperatures. We determined that our partition function is at least 95% converged at 3000 K and much better than 
this at lower temperatures. Therefore te mperatures up to 3000 K were considered. Partition function values for the parent 
isotopologue from CDMS and HITRAN ( Laraia et al.ll201ll l are compared to the present work in Table [8] 

At lower temperatures the CDMS partition function values are expected to be the most accurate and our values compare 
very well in this case. At higher temperatures only partition function values from HITRAN are available and our values are 
noticeably lower in this case. As the HITRAN partition function values are derived from using analytical approximations 
rather than a direct sum over energy levels, and do not agree as well with the CDMS values at lower temperatures, the results 
from this work are expected to be more accurate. 

For ease of use, we fitted our partition functions, Q, to a series expansion of the form used by Vidler fc Tennyson f 2000i ): 


6 

logio Q( T ) = '^2<i rl [log 10 T} n 

n=0 

with the a n values given in Table [9] 


(14) 


Table 9. Fitting parameters used to fit the partition functions. Fits are valid for temperatures between 500 and 3000 K. 



12 C 32 S 

12 C 33 S 

12 C 34 S 

12 C 3b S 

13 C“S 

13 C 33 S 

13 C 34 S 

13 C 3b S 

c»o 

-80.005915 

-76.155662 

-80.005983 

-80.004802 

-79.895945 

-79.450814 

-79.753470 

-72.151123 

ai 

135.900176 

130.515908 

135.900484 

135.901969 

136.196073 

136.926482 

136.119795 

132.718091 

CL2 

-92.069670 

-88.464147 

-92.068939 

-92.067292 

-92.226365 

-92.982691 

-92.125142 

-99.251755 

CL2, 

32.390529 

31.162671 

32.390765 

32.392690 

32.426480 

32.752272 

32.286797 

39.931597 

a4 

-6.167993 

-5.950792 

-6.168046 

-6.170901 

-6.170418 

-6.233068 

-6.095478 

-9.136353 

as 

0.599177 

0.581747 

0.599173 

0.600214 

0.599063 

0.603143 

0.582017 

1.139224 

ag 

-0.022766 

-0.022418 

-0.022766 

-0.022885 

-0.022771 

-0.022683 

-0.021361 

-0.060799 
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Table 10. Summary of our line lists. 



i2 C 32 S 

i2 C 33 S 

i2 C 34 S 

12 C 3b S 

i3 C 32 S 

13 C 33 S 

i3 C 34 S 

13 C 30 S 

Maximum v 

49 

49 

49 

49 

49 

49 

49 

49 

Maximum J 

258 

258 

258 

258 

258 

258 

258 

258 

Number of lines 

548312 

550244 

554898 

560733 

577885 

581375 

584485 

590320 


Table 11. Extract from start of states file for 12 C 32 S. The full file can be downloaded from 
http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/MNRAS//xxx/yy or www.exomol.com. 


n 

E 

9 

J 

V 

i 

0.000000 

1 

0 

0 

2 

1.634164 

3 

1 

0 

3 

4.902459 

5 

2 

0 

4 

9.804822 

7 

3 

0 

5 

16.341155 

9 

4 

0 

6 

24.511332 

11 

5 

0 


n: State counting number; 

E: State energy in cm -1 ; 
g : State degeneracy; 

J: State rotational quantum number; 
v: State vibrational quantum number. 


2.4 Line-List Calculations 

Line lists were calculated for all stable isotopologues of CS ( 12 C 32 S, 12 C 33 S, 12 C 34 S, 12 C 36 S, 13 C 32 S, 13 C 33 S, 13 C 34 S and 
13 C 36 S). These line lists span frequencies of up to 50,000 cm -1 . A summary of each line list is given in Table ITOl All rotation- 
vibration states up to v = 49 and J = 258, and all transitions between these states satisfying the dipole selection rule A J = ±1, 
were considered. 


3 RESULTS 


The line lists contain over half a million transitions eac h. For compactness an d ease of use they are separated into energy state 
and transitions files usin g the standard E xoMol format ( Tennyson et al.l2013h . which is based on a method originally developed 
for the BT2 line list I Barber et alJfioO^ L Extracts from the start of the 12 C 32 S files are given in Tables fill and fl2l The full line 
lists for all isotopologues considered can be downloaded from CDS, via ftp://cdsarc.u-strasbg.lT/pub/cats/J/MNRAS/xxx/yy, 
or http://cdsarc.u-strasbg.fr/viz-bin/qcat7J/MNRASZ/xxx/yy or can be obtained fr o m ww w.exomol.com. 

Table fl3l compares our CS line lists with the previous o nes from Pineiro et al.l (1987) and CDMS; note that the 2078 
CS transitions given in HITRAN-2012 ( Rothman et al. 2013 1 are reproduced from CDMS. Although this is an assessment 
of the quantity of the data, not its quality, it demonstrates the reason for computing the new line lists, to provide a more 
comprehensive coverage of the problem. 

Figure [2] provides an overview of CS absorption at infrared wavelengths as function of temperature. At 300 K, and below, 
the spectrum is dominated by a series of vibrational bands starting with the fundamental band at about 8 fim. At higher 
temperatures the bands become very much broader and their peak absorption is reduced. 


Table 12. Extracts from the transitions file for 12 C 32 S. The full file can be downloaded from 
http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/MNRAS//xxx/yy or www.exomol.com. 


F 

I 

A F i 

2 

1 

1.7471E-06 

3 

2 

1.6771E-05 

4 

3 

6.0640E-05 

5 

4 

1.4904E-04 

6 

5 

2.9766E-04 

7 

6 

5.2215E-04 


F: Upper state counting number; 
/: Lower state counting number; 
Api : Einstein A coefficient in s~ 1 . 
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wavelength [um] 


Figure 2. 12 C 32 S absorption spectrum at infrared wavelengths as calculated in this work for a temperature of 300 K (narrow features) 
and 3000 K (broad features). 


Table 13. Comparison of CS rotation-vibration line lists. Given are the number of isotopologues considered, the maximum values for 
vibrational (v) and rotational (J) states considered, the maximum change in vibrational state (Art) and whether intensity information 
and a partition function are provided. 


Reference 

Pineiro et al. 11987) 

CDMS 

This work 

Isotopes 

4 

6 

6 

maximum v 

20 

4 

49 

maximum J 

200 

99 

258 

maximum Av 

4 

2 

50 

Intensities? 

Yes 

Yes 

Yes 

Partition Function? 

No 

Yes 

Yes 


Comparisons with the CDMS rotational, v' — v" = 1 - 0 and v' — v‘ 
The agreement is excellent for both frequency and intensity. 


2-0 lines for 12 C 32 S are presented in Figure [3] 


Ram et al. (1995) give two figures showing their observed spectrum, a compressed view of the vibration-rotation bands 
(1000 - 1400 cm -1 ) and a portion of the R-branch region (1290 - 1310 cm~ 1 ). Emission cross- sections for 12 C 32 S were simulated 
using a Gaussian line shape profile with HWHM = 0.01 cm -1 as described in Hill et al. ( 20131 ). The resulting synthetic emission 
spectra are compared to the experimental spectra in Figures [4] and [5] For the former we the band structure and intensity 
ratio in our calculated spectra is very similar to the experiment. For the latter there is generally good agreement; however the 
intensities of five strongest lines are almost 50 % larger in the theoretical spectrum which may be due to saturation effects in 
the measured spectrum. 


4 CONCLUSIONS 

In the present work we have computed comprehensive line lists for all stable isotopologues of carbon monosulphide. We 
determined a PEC using LEVEL and modified potential parameters from the literature. We then substituted calculated 
energies in the states hie with energies derived directly from experimental frequencies to match the experimental accuracy. 
This accuracy should extend to all predicted transition frequencies up to at least v = 9 and J = 106 for 12 C 32 S, and v 
= 5 and J = 80 for 13 C 32 S, the experimental ranges. Based on comparisons with other experiments the frequencies for the 
remaining isotopologues should be predicted to sub-wavenumber accuracy at least for v < 3 and J < 21. Einstein A coefficients 
were computed from a dipole moment curve taken from the literature. Comparisons with the semi-empirical CDMS database 
suggest that the pure rotational, v' — v" = 1 - 0, and v' — v" = 2 - 0 intensities are accurate. 

The results are line lists for rotation-vibration transitions within the ground states of 12 C 32 S, 12 C 33 S, 12 C 34 S, 12 C 36 S, 
13 C 32 S, 13 C 33 S, 13 C 34 S and 13 C 36 S, which should be accurate for a range of temperatures up to at least 3000 K. The line lists 
can be downloaded from CDS or from www.exomol.com. 

Finally we note that, although our line lists are more comprehensive, for the purposes of high resolution radio astronomy 
and far-infrared studies of the low temperature objects, the CDMS line lists are recommended. 
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wavenumber [cm' 1 ] 



1200 1220 1240 1260 1280 1300 1320 1340 
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Figure 3. Absorption lines of 12 C 32 S at 300 K: ExoMol verses CDMS 
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Figure 4. 12 C 32 S at 2300 C: Emission spectrum upper, Ram et al. (1995); Emission lines lower, ExoMol. [Reprinted from Ram et al. 
(1995)] 



Figure 5. Emission spectrum of 12 C 32 S at 2300 C : upper, Ram et al. (1995); lower, ExoMol. [Reprinted from Ram et al. (1995)] 
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